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Abstract. 

We present time-dependent solutions of thin, supersonic accretion 

flows near a black hole and compare them with analytical solutions. Such 

^ ■ flows of inviscid, adiabatic gas are characterized by the specific angular 

'sj" \ momentum and the specific energy. We confirm that for a wide range of 

|0 ■ above parameters a stable standing shock wave with a vortex inside it 

forms close to the black hole. Apart from steady state solutions, we show 
(^ ■ the existence of non-steady solutions for thin accretion flows where the 

[~^ i accretion shock is destroyed and re-generated periodically. The unstable 

0^ ■ behavior should be caused by dynamically induced instabilities, since 

inviscid, adiabatic gas is considered. We discuss possible relevance of 
Qh' the periodic behavior on quasi-periodic oscillations (QPOs) observed in 

Q . galactic and extragalactic black hole candidates. 



Introduction 



Rotating accretion flows are important ingredients in many astrophysical sys- 
C^ I tems containing a black hole, which involve mass transfer from one object to 

another (such as in a binary system) or from set of objects to another (such 
as in a galactic center). The standard disk model of such accretion flows by 
Shakura & Sunyaev (1973) assumes Keplerian distribution of accreting matter. 
There, the inner edge of the disk is chosen to coincide with the marginally stable 
orbit located at three Schwarzschild radii, r, = 2>Rg, where the Schwarzschild 
radius, Rg = 2GMbh /(? ■, is the horizon of a black hole of mass Mbh- This disk 
model is clearly incomplete, since the inner boundary condition on the horizon 
was not taken care of. As an accretion flow approaches the horizon, its radial 
velocity reaches the velocity of light. Therefore, a black hole accretion flow is 
necessarily supersonic and must pass through a sonic point where the flow has 
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to be sub-Keplerian. Thus, independent of heating and coohng processes, a 
black hole accretion has to deviate from a standard Shakura-Sunyaev type Ke- 
plerian disk. The disk with realistic accretion flows is called the advective disk 
(Chakrabarti 1996). Here, we present the results of numerical study of accretion 
flows near a black hole by assuming they are thin, axisymmetric, and inviscid. 

2. Analytic Consideration 

We choose cylindrical coordinates (r, 9, z) and place a black hole at the center. 
We assume that the gravitational field of the black hole can be described in 
terms of the potential introduced by Paczyhski & Wiita (1980) 

At \ GMhh , . 

nr,z) = -- — (1) 

It itg 

where R = \/r'^ + z^. The accreting matter is assumed to be adiabatic gas with- 
out cooling and dissipation and described with a polytropic equation of state, 
P = Kp"', where 7 is the adiabatic index which is considered to be constant with 
4/3 throughout the flow. K is related to the specific entropy of the flow, s, and 
varies only at shocks, if present. Since we consider only weak viscosity limit, 
the specific angular momentum of the accretion flow, A = vvq, is assumed to be 
conserved. Thus, unlike a Bondi flow, which is described by a single parameter 
(say, specific energy), the one-dimensional accretion flows are described by two 
parameters which are the specific energy, £, and the specific angular momen- 
tum, A. Fig. 1 shows the classification in the parameter space (Chakrabarti 
1989; Ryu & Chakrabarti 1997): 

N: No sonic points. Shock only if supersonic injection. 
O: Outer sonic point only as in a Bondi solution. No shock. 
I: Inner sonic point only. Shock only if supersonic injection. 
O*: Outer and center sonic points. Solution does not extend to the horizon. 
I*: Inner and center sonic points. Solution does not extend to large distance. 
SA: Two (outer and inner) sonic points. Shock in accretion solutions but not in 
wind solutions. 

SW: Two sonic points. Shock in wind solutions but not in accretion solutions. 
NSA: Two sonic points. No shock condition satisfied in accretion solutions. 
NSW: Two sonic points. No shock condition satisfied in wind solutions. 

Although these solutions are strictly valid for inviscid flows, even when 
viscosity is high, given that the viscous time-scale is likely to be much larger 
compared to the infall time-scale, the inviscid solutions are likely to remain 
important. 

3. One-dimensional Numerical Solutions 

Fig. 2 shows an example numerical solution from the 'SA' region. We superposed 
analytical solution (solid lines) and numerical solutions with TVD code (dashed 
lines) and with SPH code (dotted lines) (Molteni, Ryu & Chakrabarti 1996). 
The TVD calculation was done with 512 grids and the SPH calculation was 
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Figure 1. Classification of one-dimensional accretion flows in the pa- 
rameter space of specific energy, iS, and the specific angular momentum, 
A. E is given in unit of (? and A in unit of RqC. 



1.5Rg for 
1.25Rg for the SPH calculation. The adiabatic 



done with ~ 560 particles of size h = 0.3Rg. Matter was injected at the outer 
boundary located at x = 50Rg. An absorption condition was used to mimic the 
black hole horizon at the inner boundary which is located either at x 
the TVD calculation or at x 
index 7 = 4/3 was used. 

The figure shows an excellent agreement between the analytic and numerical 
solutions. Here, the flow starts out subsonically, presumably from a Keplerian 
disk. Then, it enters through the outer sonic point (located at x = 27.97Rg), 
passes through the shock (located at x = 7.98Rg), and subsequently passes 



through the inner sonic point (located at x 
hole. 



2.57Rg) before entering the black 



4. Two-dimensional Numerical Solutions 

In multi-dimensional accretion flows, non-steady solutions, as well as steady 
solutions, exist. Here, we discuss some examples of accretion flows with zero 
specific energy {£ = 0). More extensive discussion on accretion flows with £ = 
was reported in Ryu, Chakrabarti & Molteni (1997). Discussion on extensive 
calculations for accretion flows with non-zero specific energy {£ 7^ 0) will be 
reported elsewhere (Ryu & Chakrabarti 1997). 

In these simulations, we inject supersonic matter (with a radial Mach num- 
ber M = Vr/a = 10) at the outer boundary, rf, = 50Rg. The inflow at the 
outer boundary is assumed to have a small thickness, hm, or a small arc angle. 
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Figure 2. One-dimensional accretion flow with a standing shock in 
the region of 'SA'. £ = 0.036 and A = 1.80 are used. The solid line 
is the analytical solution, and the long and short dashed lines are the 
solutions of the TVD and SPH simulations, respectively. Upper panel 
is the mass density in arbitrary units and the lower panel is the Mach 
number of the flow. 

9in = aicctan{hin/rh) <^ 1. If zero-energy accretion flows belong the region of 
'O', most of the material is accreted into the black hole forming a stable quasi- 
spherical flow or a simple disk-like structure around it (just as a Bondi flow). If 
accretion flows belong the region of 'SA' and 'O*', the incoming material pro- 
duces a stable standing shock with one or more vortices behind it and some of 
it is deflected away at the shock as a conical outgoing wind of higher entropy. 

Accretion flows with parameters in the region 'NSA' with 1.782RgC < A < 
1.854:RgC show an unstable behavior. Fig. 3 shows an example numerical so- 
lution with A = 1.85RgC. In this case, the structure with an accretion shock 
and a generally subsonic high density disk is established around the black hole. 
However, the structure is not stable. At the accretion shock, the incoming flow 
is deflected. But some of the post-shock flow, which is further accelerated by 
the pressure gradient behind the shock, goes through a second shock, where the 
flow is deflected once more downwards. The downward flow squeezes the in- 
coming material, and the accretion shock starts collapsing (t = 1.4 x W^Rg/c). 
In the process of the collapse, some of the post-shock material escapes as wind 
but most is absorbed into the black hole. After the collapse, the re-building of 
the accretion shock starts with the incoming material bouncing back from the 
centrifugal barrier. The subsonic post-shock region becomes a reservoir of ma- 
terial, so the material is accumulated behind the shock. With the accumulated 
material a giant vortex is formed, which in turn supports the accretion shock 
(i = 2 X W^Rg/c). This continues until the incoming flow is squeezed enough so 
the accretion shock collapsed, and the cycle continues. 
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Figure 3. Two-dimensional accretion flow showing an unstable be- 
havior in the region of 'NSA'. £ = and A = 1.85 are used. Contours 
are for density and arrows are for velocity vector. Time is given in unit 
of Rg/c. 



5. Discussion 

The time scale of the periodicity of unstable flows is interesting. It is in the 
range of 

^P,4-6xl03:^ = 4-6x 10-2 (^] s. (2) 

The modulation of amplitude is also very significant, and could be as much as a 
hundred percent depending on detailed processes. Oscillations with these char- 
acteristics have been observed in black hole candidates and are called the QPOs. 
For instance, in the QPOs from the low mass x-ray binaries, the oscillation fre- 
quency has been found to lie typically between 5 and 60 Hz (Van der Klis 1989). 
Thus, compact objects with mass M ~ 0.3 — 5Mq could generate oscillations 
of right frequencies due to the instability discussed in this paper. Similar os- 
cillations of period on the order of a few hours to a few days are expected in 
soft X-rays and UV emissions from galactic centers. However, by considering 
simplified physics which we have assumed here, it may be premature to assume 
that the presented mechanism would explain all the QPOs observed. In some 
of the cases, the oscillation may be due to the dynamic instability considered 
in Ryu et al. (1995), or it may be due to resonance of the cooling time scale 
(bremsstrahlung or Comptonization, whatever the case may be) and the infall 
time scale in the enhanced density region near the centrifugal barrier as shown 
by Molteni, Sponholz and Chakrabarti (1996). In detailed works, radiative pro- 
cesses as well as viscosity should be included in the accretion calculations to 
examine the observational consequences of the present instability. 
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